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ABSTRACT 

We show that a < 0, (3 > 0,7 = 5 = Painleve III equation arises as a zero- 
curvature condition in the Belinskii-Zakharov inverse scattering formulation for 
Bianchi cosmological models. For special values of the parameters this Painleve 
III equation becomes the dynamical equation for Bianchi I, II, VIq and VIIq models. 
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1. Introduction 



Belinskii and Francaviglia showed in [1] that the Einstein equations for Bianchi 
I, II, VIo and VIIo spacetimes admit a zero-curvature representation, i.e. they 
found a linear system whose integrability condition is the dynamical equation 
for the Bianchi model. This was done by using a more general framework of 
Belinskii-Zakharov (BZ) inverse scattering method for the spacetimes admitting 
two commuting spacelike Killing vectors [2] . The results of [1] demonstrated that 
the Bianchi models which admit two commuting spacelike Killing vectors are solv- 
able dynamical systems. However, not much work has been done on the issue of 
what kind of integrable nonlinear dynamical equations can be obtained from this 
approach. 

In [3] it has been shown that in the case of the Bianchi VIIo model one obtains 
a special Painleve III (PHI) equation, which is a = — 2, (3 = 2, ^ = 5 = case of 
the standard PHI form [4] 



In this paper we consider the Belinskii-Zakharov formulation for all Bianchi models 
which admit two commuting Killing vectors, i.e. types I, II, VIo and VIIo- We 
only consider the equations of motion for local degrees of freedom, and we do not 
discuss the problems related with non-trivial topology of the spatial manifold [3]. 
By considering a larger class of Bianchi models, we obtain a more general PHI 
equation as the dynamical equation, namely a < 0, (3 > 0, 7 = 5 = 0. Although 
the relevant Bianchi models correspond to special values of the parameters, we 
show that the zero curvature representation is valid for all other values of the 
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parameters. Consequently we obtain a zero- curvature representation of a < 0, 
/3 > 0, 7 = 5 = Painleve III equation in the Belinskii-Zakharov inverse scattering 
formulation. 



2. Zero-curvature formulation for Bianchi models 

Belinskii-Francaviglia approach to solving the dynamics of Bianchi models [1] 
is derived from the Belinskii-Zakharov method for solving the Einstein equations 
for spacetimes with two commuting spacelike Killing vectors [2]. Such spacetimes 
have the following form of the metric 



where a,b — 1,2, {x°, x 1 , x 2 , x 3 } = {t, x, y, z}, f is a positive function and g a \, is a 
symmetric two-by-two matrix. It is convenient to introduce the null co-ordinates 
(£->v) — {\{ z + *)) \{ z ~~ t)), since the form of the metric (2.1) is preserved by the 
conformal co-ordinate transformations (£,rj) — > (£(£), viv))- The positivity of the 
function / is preserved if d v fj > 0. 

The complete set of vacuum Einstein equations for the metric (2.1) decomposes 
into two groups of equations [2]. The first group determines the matrix and 
can be written as a single matrix equation, called the Ernst equation 



where a 2 = detg. The second group of equations determines the function f(£,r)) 
in terms of a given solution of the Ernst equation 



ds 2 = f{t, z){- dt 2 + dz 2 ) + g ab (t, z) dx a dx b , 
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dl(\na) 1 9 
In/ = + -. tr5 2 , 2.4 

where = d^a, a n = d v a and the matrices A and B are defined by 

A — —a % g~\ B = a d n g g~ l . (2.5) 

Thus the dynamics of the system is determined by the Ernst equation (2.2). An 
important consequence of the Ernst equation is that = so that a = c(£)+d(r)). 
By using the conformal transformations, one can bring the functions c(£) and d(rj) 
to a prescribed form. 

The crucial step in the inverse scattering method is to define the linearized 
system whose integrability conditions are the equations of interest, in our case the 
equation (2.2). Following the ref. [2], we define two differential operators 

D 1 = dt -^±d x ,D 2 = d v+ ^d x , (2.6) 
A — a A + a 

where A is a complex parameter independent of the co-ordinates {£,77}. The dif- 
ferential operators D\ and D2 commute since a satisfies the wave equation, and 
hence one can consider the following linear system 

Dnl> = -^—il> , D 2 ^ = -^—iIj , (2.7) 
a — <J A + a 

where £,77) is a complex matrix function. The integrability condition for the 
system (2.7) is given by the Ernst equation (2.2). Furthermore, a solution *0(A, £, 77) 
yields a matrix g(^ : rj) that satisfies the Ernst equation (2.2). Namely, the matrix 
g(t,v) is given by 
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(2.8) 



In order to take into account that g(£, rj) is real and symmetric we have to impose 
two additional conditions, see [2]. Also, it is easy to see that the equations (2.7) 
for A = 0, imply equations (2.5). 

Bianchi spacetimes (see [5] for a review and references) have finitely many 
degrees of freedom, and only the Bianchi types I, II, VIo and VIIo admit two 
commuting spacelike Killing vectors. The metric for these Bianchi spacetimes has 
the form (2.1), and this can be shown by considering the general Bianchi spacetime 
metric 

ds 2 = -dT 2 + gij (T, x k ) dx i dx j , (2.9) 

where gij = gij{T) x 1 ii xk ) X J j( xk ) an d X T ( xk ) are the one-forms associated with 
the spatial manifold. These one-forms satisfy the Maurer-Cartan equations 

dx 1 + h cI JKX J ^X K = , 

where the structure constants C 1 jk correspond to the Lie algebra of the symmetry 
group of the Bianchi model. For the relevant models the structure constants satisfy 

C I j K = ej KL S LI , (2.10) 

where cjkl is a totally antisymmetric tensor density and S is a symmetric matrix. 
In this case the one forms x 1 take the following form 

X 1 = l\(z) dx + l\(z) dy , X 2 = dx + l 2 2 (z) dy , X 3 = dz .(2.11) 

An important consequence of the Maurer-Cartan equations for the one forms x 1 is 
that the matrix I = \ \l a b\ \ satisfies the following linear differential equation 



where the matrix C is the upper two-by-two block on the principal diagonal of the 
matrix S IJ and e is the antisymmetric matrix with ei2 = 1. 

After a time redefinition t = t(T), the metric (2.9) can be written in the form 
(2.1) 

ds 2 = fit) (-dt 2 + dz 2 ) + g ab (t, z) dx a dx b , (2.13) 

where 

g(t,z) = l T (zmt)l(z) , (2-14) 

and 7 is a two-by-two symmetric matrix. Notice that now a 2 = (det I) 2 det 7, and 
since det / = 1 we get 

a 2 (t) = det 7(f) . (2.15) 

In addition, a has to satisfy the wave equation, so that a(t) = 0, and hence a can 
only be a linear function of time. 

The linearized system (2.7) can be simplified for the models described by the 
metric (2.13). The first step is to define a two-by-two matrix function ip by 

ip(t, z, A) = l T {z) ip(t, z, A) l(z) , (2.16) 

and a constant two- by-two matrix R — e C. The second step is to substitute (2.14) 
into (2.5) and use the definition of the coordinates £ and rj. Then the results of 
these calculations, together with the definition (2.16), can be used to simplify the 
equations (2.7). The crucial step in which a simplification occurs is to perform 
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a conformal coordinate transformation {t, z, A} — > {t, u>,A}, where w is given by 
u> = \{\ + 2(3 + A). The linear system after this coordinate transformation 
involves only derivatives in t and A since all the terms involving derivatives in w 
are canceled. Finally, it is useful to make some simple linear combinations of the 
two equations and to use the fact that a is a linear function of time. In this way 
one obtains a new linear system 

dt<p = j {iR T r\ - vR T ) , 

d x p = X - (-Rip - pR T + ^77" V + ^R T - ^lR T r\) , (2.17) 

where we have set a = t. 

Although the matrix function ip(t,X,w) depends on all three variables, the 
right-hand side of the system (2.17) does not have any w dependence. The inte- 
grability condition for the system (2.17) is 

-^-(i77 _1 ) =RiR T T X -lR T l~ 1 R ■ (2-18) 

Equivalently, one can derive the equation (2.18) by a direct substitution of the 
formula (2.14) into equation (2.2). Thus the dynamics of these Bianchi models is 
determined by the equation (2.18). 
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3. Zero-curvature representation for Painleve III 



The linear system (2.17) and the corresponding nonlinear equation (2.18) were 
derived for special matrices R, which correspond to C matrices of the relevant 
Bianchi models. Note that the C matrix is symmetric, and the only relevant 
information about the Bianchi model is contained in its signature, so that there 
are four distinct possibilities: 

(1) C — diag(l, 1) for Bianchi VIIo, 

(2) C = diag(l, -1) for Bianchi VI , 

(3) C = diag(l, 0) for Bianchi II, 

(4) C = diag{0, 0) for Bianchi I. 

However, if we consider the linear system (2.17) independently of Bianchi models, 
then we can take C to be an arbitrary symmetric two-by-two matrix. In this case 
we have 




(3.1) 



with R = eC, 7 = diag(a, b) and ab = t 2 



Then the consistency condition (2.18) 



gives 



t 1 — (ta 1 d) = k 2 a 1 b — c 2 ab 



dV ' 

= d{k + cab- 1 ) . 



-i 



(3.2) 



The second equation implies either d = or a 2 = —{k/c)t 2 . The second possibility 
gives a linear in time solution, and the first possibility is more interesting, d — 
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case gives 

t _1 ^(fa _1 d) = k 2 a- 2 t 2 - c 2 a 2 t- 2 . (3.3) 

By making a change of variables u = a 2 /t 2 , r = t 2 /A, the equation (3.3) takes the 
standard Painleve III form (1.1) 

d 2 u 1 fdu\ 2 1 fdu\ 2/ 9 9 , 9 \ . 

- 7 y+ 7 (-^ 2 +* 2 )' (3.4) 

so that a = -2c 2 , (3 = 2k 2 and 7 = 5 = 0. 

Let us now consider the Bianchi models. In the case of Bianchi VIIo model 
the spatial hyper-surface is a three torus T 3 . The matrix C is diag{l, 1) and the 
matrix / is given by 

(cos z sin z \ 
(3.5) 
— sin z cos z J 

The matrix R is given by 

fl =(-°i o) • (X6) 

so that the local dynamics is given by (3.3) for c = k — 1. By making a change of 
variables u = e q , this PHI equation takes a more symmetric form 

4-(t^] =-4sinhg . (3.7) 
dr \ dr J 

Bianchi VIo model corresponds to C = diag(l, —1), so that the equation (2.12) 
gives 

(cosh z sinh z \ 
sinnz cosh z / 

This model has a non-compact spatial manifold, which is locally compact, and the 
local dynamics is the same as in the Bianchi VIIo case, because c 2 = k 2 = 1. 



In the Bianchi II case C = diag(l, 0), and the dynamics is given by (3.3) with 
c = 1 and k — 0. The I matrix is given by 

K*)=(^ Q , (3-9) 

and this model allows compact spatial sections [6]. In the Bianchi I case C = 0, so 
that c = k = and I = diag(l, 1). 



4. Conclusions 

We have shown that a class of PHI equations (a>0,/3<0,7 = 5 = 0) arises 
as a zero-curvature condition in the Belinskii-Zakharov inverse scattering method 
applied to Bianchi spacetime metrics. For the particular values of the parameters 
these PHI equations become the dynamical equations for Bianchi I, II, VIo and VIIo 
models. Note that the linear system (2.17) can be transformed into the standard 
form 



ox= A * ■ m =m (41) 



where A and B are four-by-four matrices given by 



A=\ (7i? T 7 _1 ® h + h <8> R T ) 
A 

B = \{-R+ jil- 1 - ^7^ T 7 _1 ) ® h + \ - lj h ® R , (4.2) 



where I2 is the identity matrix and \P is a column formed from the columns of the 
matrix </?. 
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This linear system is different from the linear system which is used for the study 
of Painleve III equation within the isomonodromic deformation (IMD) method 
[7,8,9]. Although it is not obvious what are the advantages of the new linear system 
in comparison with the IMD linear system, there are some interesting features of 
the new system which can be investigated. 

In the Belinskii-Zakharov inverse scattering approach, it is natural to consider 
path-ordered exponentials of matrices (holonomies) associated to the linear system 
(2.7) in order to find the integrals of motion [10], [11]. In the special case of Bianchi 
metrics, the holonomy construction would simplify, and one could try to see what 
kind of expressions one would obtain for the PHI equation. 

Note that the new linear system can be interpreted as a Lax pair for a dynam- 
ical system with a time- dependent Hamiltonian 

H = —p 2 + 2(k 2 e- q + cV) . (4.3) 
2r 

This is analogous and complementary to the results of Harnad and Routhier [12], 
where a Lax pair for PHI equation with afijS ^ was constructed. It is interesting 
that in that case the Lax pair contains f dm, which does not happen in our case 
of PHI equation. 

In the context of Bianchi models, it is more natural to work with a dynamically 
equivalent Hamiltonian to (4.3) 

H = \p 2 + 2e i {k 2 e- q + cV) , (4.4) 

where r = e t . One can now examine the physical properties of the solutions, like 
small and large time asymptotic, as well as the singularities, since these properties 
of the Painleve III solutions have been thoroughly studied [9] . 
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